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Abstract: We present an updated dynamical and statistical analysis of outer Oort cloud 
cometary evidence suggesting the sun has a wide-binary Jovian mass companion. The results 
support a conjecture that there exists a companion of mass « 1 — 4 MjupUer orbiting in the 
innermost region of the outer Oort cloud. Our most restrictive prediction is that the orientation 
angles of the orbit normal in galactic coordinates are centered on fl, the galactic longitude of 
the ascending node = 319° and i, the galactic inclination = 103° (or the opposite direction) 
with an uncertainty in the normal direction subtending Ri 2% of the sky. A Bayesian statistical 
analysis suggests that the probability of the companion hypothesis is comparable to or greater 
than the probability of the null hypothesis of a statistical fluke. Such a companion could also 
have produced the detached Kuiper Belt object Sedna. The putative companion could be easily 
detected by the recently launched Wide- field Infrared Survey Explorer (WISE). 

Key Words: Comets, dynamics; Celestial mechanics; Kuiper Belt; Jovian planets; Planetary dynamics 
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1. Introduction 



Anomalies in the aphelia distribution and orbital elements of Outer Oort cloud comets led to the 
suggestion that fti 20% of these comets were made discernable due to a weak impulse from a bound Jovian 



mass body (jMatese et al 



19991 )). Since that time the data base of comets has doubled. Further motivation 



for an updated analysis comes from the recent launch of the Wide-field Infrared Survey Explorer (WISE; 



Wright 



2007))), which could easily detect the putative companion orbiting in the outer Oort cloud. Such 



an object woul d be incapable of creating com e t "storms" . To help mitigate popular confusion with the 



Nemesis model ( Whitmire and Jackson 



1984 . 



Davis et al. 



1984 )) we use the name recently suggested by 



Kirkpatrick and Wright (2010), Tyche, (the good sister of Nemesis) for the putative companion. 

The outer Port cloud (OOC ) is formally defined as the ensemble of comets having original semimajor 



axes A > 10^ AU (Oort 



19501 )). It has been shown that the majorit y of these c omets that are 



made discernable are first-time entrants into the inner planetary region ([Fernandez 



19811) 1 and these 



comets are therefore commonly referred to as new. The dominance of the galactic tide in making OOC 



comets discernable at the present epoch has been predicted on theoretical grounds (IHeisler and Tremaine 



19861)'). Observational evide n ce of this dominance has been claimed to be compell i ng 



Matese and Whitman 



1993); 



Wiegert and Tremaine 



19991) : 



Matese and Lissauer 



Delsemme 



20041) ') 



IMS); 



Matese and Lissauer (2004) adopted an in situ energy distribution similar to the initial distribution of 
Rickman et al. (2008) and took the remaining phase space external to a "loss cylinder" to be uniformly 
populated at the present epoch. The distribution of cometary orbital elements made discernable from the 
tide alone was then obtained and compared with observations. 



Similar modeling (jMatese and Lissauer 



20021 )) had been performed including single stellar impulses 
which mapped the comet flux over a time interval of 5 Myr, in 0.1 Myr intervals. Peak impulsive 
enhancements > 20% were found to have a half-maximum duration of « 2 Myr and occurred with a mean 



time interval of ~ 15 Myr. Various time- varying distributions of elements were compared with the modeled 
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tide-alone results and inferences about the signatures of a weak stellar impulse were drawn. 



In Section 2 we review a discussion ( Matese and Lissauer 



20041) ) of a subtle characteristic of galactic 



tidal dominance which is difhcult to mimic with observational selection effects or bad data. Along with 
the more well known feature of the deficiency of major axis orientations in the direction of the galactic 
poles and equator, we compare with observations these predictions based on the tidal interaction alone and 
show that the data are of sufficiently high quality to unambiguously demonstrate the dominance of the 
galactic tide in making comets discernable at the present epoch. A critiq ue of objections to this a ssertion 



(Rickman et al. 



20091) ) provide 



20081 )) is also presented. More recent detailed modeling (|Kaib and QuinnI ( 
important insights into the evolving populations of the in situ and discernable populations of the Oort 
cloud. We comment further on these works in this section. 

In Section 3 we describe the theoretical analysis combining a secular approximation for the galactic tide 
and for a point mass perturber, describing how a weak perturbation of OOC comets would manifest itself 
observationally. Evidence suggesting that ther e is such an al i gned i mpulsive component of the observed 



OOC comet flux has been previously reported (jMatese et al. 



known observational biases can explain the alignment found there (jHorner and Evans 



1999 )). It has been found tha t none of the 



20021 )). The size of 



the available data has since doubled which leads us to review the arguments here. 

In Section 4 we present the supportive evidence that an impulsive enhancement in the new comet flux 
of w 20% persists in the updated data. We also discuss dynamical and observational limits on parameters 
describing the putative companion. Section 5 summarizes our results and presents our conclusions. 



2. Secular dynamics of the galactic tide 



Near-parabolic comets are most likely to have their perihelia reduced to the discernable region. The 
dynamics o f the galactic t i de acti n g on near-parabolic OOC co m ets is most simply described in a Newtonian 



framework 



Matese et al. 



1 



19991) ■ iMatese and Lissauer 



20021) . 



Matese and Lissauer 



20041 )). A summary 
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of their analyses is now given and followed with the evidence that the galactic tidal perturbation dominates 
in making OOC comets discernable at the present epoch. 



2.1. Theory 

Saturn and Jupiter provide an effective dynamical barrier to the migration of OOC comet perihelia. 
OOC comets that are approaching the planetary zone at the present time were unlikely to have had a 
prior perihelion, Qprior, that was interior to the "loss cylinder" radius, qic ~ 15 AU, when it left the 
planetary rGgion on the present orbit. The simplifying cLSSumption qiq ^ Qprior 

is then made for the 

present orbit. During the present orbit comet perihelion will then have been changed by the galactic tide 
(and by any putative companion or stellar perturbation). The orbital elements just before re-entering 
the planetary region on the present orbit are commonly referred to as "original" and will be, in essence, 
the observed values with the exception of the semimajor axis (perturbations by the major planets do not 
significantly change any other orbital element of OOC comets). Thus qprior is changed to Qoriginai ~ lobs, 
the observed value, during the course of the present orbit. As an observed comet comes within a discernable 
region (gobs < Qdiscemabie ~ 5AU) and leaves the planetary region again, the semimajor axis will have 
been changed from Aoriginai to Ajuture, i-e., Afuture = Aprior for the next orbit. The comet is most 
likely ejected or turned into an inner Oort cloud comet with a small fraction returning as OOC comets. 
Daughter comets return i ng to the discernable zone are likely to have faded and be more difficult to observ e 



( Wiegert and Tremaine 



20081) ') 



19991) ). The Ylth Catalogue of Cometary Orbits ([Marsden and Williamsl 
indicates that 14% of observed original OOC comets exit the planetary region as future OOC comets. 
Therefore the discernable population of OOC comets should be dominated by first time entrants to the loss 
cylinder. In the following we adopt the notation A = Aoriginai- 

For near-parabolic comets, the angular momentum per unit mass determines the perihelion distance. 



H = R X R, (H _L q, iJ « y/2fj,Qq). With these assumptions, an observed OOC comet entering the loss 
cylinder region for the first time had perihelion distances qobs < qdiscemabie < qic < qprior- Therefore, 
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reducing g in a single orbit requires a decrease in angular momentum from the galactic tidal torque (and/or 
from angular momentum changes by the putative companion or star), 



or, Qprtor " Qobs = (AH^ - 2Hobs ' AH) /2H(. 



(1) 



The weakest perturbation that could make a comet discernable would reduce the prior perihelion distance 

from Qpnor ~ Qlc tO Qobs ~ Qdiscernable (sCC Fig. 1) SUCh that 



\^^\rmn = \/2Ai0 Qlc ' \/2/i0 Qd: 



iscemable 



(2) 



Also of interest is the evolution of the aphelion orientation, Q = (cos B cos L, cos B sin L, sin B), expressed 
in terms of the aphelion latitude, B, and longitude, L. 

If the galactic tide dominates in making OOC comets discernable we can recast Eq. (1) as 

Qprior Qobs 

= (AHtide'^ — 2Hobs • AHtide) /2/i0 from which it is implied that if weak tidal perturbations 
dominate in making OOC comets discernable, the tidal characteristic, S = Sig nfHnhs • AHtide) will more 



often be —1 than +1. Detailed modeling results (jMatese and Lissauer 



20041 )) confirm this implication. 



This characteristic combination of observed orbital elements forms an essential aspect of the pr e sent 



analysis and has not be en included in modeling results presented elsewhere {e.g., 



Rickman et al 



Kaib and Quinn 



20091)). 



200M) 



A graphical illustration of this theme in Fig. 1 shows the phase space changes in H for a specific 
choice of Q and A over the course of a single orbit period. As comets recede from the planetary region on 
their prior orbits, the interior of the loss cylinder phase space region is essentially emptied of OOC comets 
by planetary perturbations. The adjacent exterior region remains uniformly populated in this model. To 
lowest order in q/A, i.e., the near-parabolic phase space region just outside the loss cylinder boundary, the 
vector displacement in specific angular momentum, AHtide, is independent of the prior value of angular 
momentum, Hprior, and depends only on A and the major axis orientation, Q, which are taken to be fixed 
for the phase space of Fig. 1. In a single orbit all nearby specific angu lar momentum phase space points, 



filled and empty, are displaced uniformly (jMatese and Lissauer 



20041 )). Figure 1 illustrates that in the 




Fig. 1. — Schematic illustration of the plane of specific angular momentum phase space, H, for a given A 
and Q, the aphelion direction defining the normal to this plane. Phase space points illustrated have the 
same values of A, B, and L, but cover the entire region of H just outside the loss cylinder. Light gray 
shaded region outside the loss cylinder boundary {H > Hic oc denotes the filled phase space of comets 
leaving the planetary region on their prior orbits. The uniform displacement of all phase space points one 
orbit later by a tidal perturbation AHtide is shown partially refilling the loss cylinder region (indicated 
as the dark gray region). The discernable zone within the loss cylinder {H < Hdiscemable oc -s/qdiscemable) 
is indicated. It remains unfilled for negligible perturbations of small-A comets, becomes totally refilled 
for strong perturbations of large-A OOC comets, and, as illustrated here, is partially refilled for weak 
perturbations of intermediate-^ OOC comets. The tidal characteristic S = Sign(Hobs ■ AHtide) = —1 for 
all observed comets in the refilled portion of the discernable phase space here. Over successive orbits the 
loss cylinder and its shadow region of unfilled phase space (shaded white) will be depleted for small-A and 
intermediate-A OOC comets. 



present modeling, the phrase "loss cylinder" might be more appropriately changed to "loss circle" . It also 
indicates that the phrase "loss cone" which is commonly used in describing impulsive domination is no 
longer appropriate if the tide dominates. 

The magnitude of the single-orbit angular momentum displacement is strongly dependent on A, varying 
as A''/^. Small-A comets are defined here as having negligible galactic perturbations, < |AH|^.^, unable 
to repopulate the discernable zone with new comets. This includes the inner Oort cloud (IOC), but also 
overlaps with the smallest- A population of the formal OOC. Large-^ OOC comets are defined here as those 
having strong tidal perturbations resulting in large displacements in H that completely refill the discernable 
zone, making S = ^1 equally likely. Intermediate- A OOC comets are defined here as those that are weakly 
perturbed by the tide and only partially refill the discernable zone. Intermediate-^ comets preferentially 
have S = —1, as seen in Fig. 1. In this context intermediate- A OOC comets have the smallest observed 
values of A among the OOC comets if the tide alone makes the comets discernable. The fuzzy boundaries 
between small-^, intermediate-^ and large-A depend weakly on B and L. We choose the boundaries for 
these intervals based on data discussed in Section 3. 

Therefore, independent of stellar influences on the in situ distribution of semimajor axes, if the galactic 
tidal interaction with the OOC dominates impulsive interactions in making OOC comets discernable at the 
present epoch we should see 

• a preponderance of OOC comets with S = —1 over those with S = +1, and 

• an association in which S = —1 correlates with the smallest observed values of A for OOC comets. 

Conversely, if perturber impulses dominate in making OOC comets discernable at the present epoch, the 
unique tidal characteristic S should be a random variable and should be uncorrelated with A. 
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2.2. Observational evidence for tidal dominance at the present epoch 

Our data are taken from the 17th Catalogue of Cometary Orbits from which we convert the echptic 
Eulerian orbital angles into the galactic angles B, L and a, the orientation angle of H defined in Matese and 
Lissauer (2004). The Catalogue lists 102 OOC comets (see Appendix) of the highest quality class, lA (we 
count the split comet C/1996-J1A(B) as a single comet ). The quality c lass pr edominantly distinguishes the 



accuracy of the original semimajor axis determination (jMarsden et al 



19781) ). Since our analysis depends 



sensitively on an accurate determinatio n of A, we restrict our detailed discussions to class lA comets. In a 



previous analysis ([Matese et al. 



19991 )) the orbital elements of 82 OOC comets of quality classes lA+lB 



given in the 11th Catalogue were used, 47 of which were class lA. 
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Fig. 2. — The distribution in | sini?| for class lA new comets. 



The first observational indi cation that the ga lactic tide do minates involved the distribution in the 



galactic latitude of aphelion, B (jPelsemme 



19991 )) that the dominant 



1987n V One finds (jMatese et al. 
disk tidal term in the perturbation is cx | sin B cos -B| . (Rickman et al. (2008) conclude that this perturbation 
is cx |sin_B| rather than cx | sin_Bcos_B|). B and L are sensibly constant in the course of a single orbit 
since Q has significant inertia for near-parabolic orbits. Therefore, if the tide dominates we sliould expect 
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deficiencies of major axis orientations along the galactic poles and tfie galactic equator, and peaks near 
B = ±45°. One migfit argue that observational selection effects can artificially produce an equatorial gap, 
but polar gaps will be more difficult to attribute to an observational selection effect. 

In Fig. 2 we show the results presented as a distribution in | sin_B|, which would be uniform for a 
random distribution. Polar and equatorial gaps are clear, as predicted if the galactic tide dominates. A 
small (but potentially informat i ve) dis crepancy is the location of the peak. If the tide dominates, our 



modeling (IMatese and Lissauer 



_2004|)) predicts a peak at | sin_B| w 0.7, somewhat larger than seen in the 

data. We now look to the tidal characteristic S to further emphasize tidal dominance. 
70 




Fig. 3. — The cumulative binding energy distribution {x = 10^ A\J / A) separately illustrated for S = ^1. 
Solid line o 5 = -1, dotted line o S* = +1. 



The prediction that S* = ^1 should be equally likely for large- A OOC comets, and that there should 
be a preponderance of S* = — 1 for intermediate-A OOC comets (see Fig. 1) is now considered. In terms of 
the original orbital binding energy parameter x = 10^ A\J /A, class lA OOC comets have a mean formal 
error of fti 5 un it s, but the un c ertain ty due to unmodeled outgassing effects is likely to be somewhat larger 



(Kresak 



1992I) . 



Krolikowska 



2006|)). In particular the number of nominally unbound original orbits is 
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likely indicative of the true errors embodied in the tails of the class lA OOC energy distribution. 

In Fig. 3 we show the cumulative class lA binding energy distribution of 66 comets with S = —1 and 
36 comets with S" = +1. The binomial probability that as many or more would exhibit this imbalance if in 
fact S = =f1 were equally likely is 2 x 10^^. Further, as predicted by tidal dynamics, this preponderance 
of S* = — 1 also correlates with intermediate- OOC comets in a statistically significant manner. S = =Fl 
is approximately equally likely for comets with a; < 30, suggesting that for this range of semimajor 
axes the galactic tide is s trong enough to refill the di scernable zone almost completely, i. e., the tidal 



efficiency is nearly 100% (jMatese and Lissauer 



20021) ). We therefore take this as the defining boundary 



between large-^ and intermediate-A comets. In contrast, for 30 < a; the evidence is that the tide weakens 
dramatically since 45 comets have S — —1 and only 12 have = +1. The binomial probability that as many 
or more would exhibit this imbalance if in fact S — =Fl were equally likely for this energy range is 7 x 10^^. 
This is unambiguous evidence that the class lA data are of sufficiently high quality and sufficiently free 
of observational selection effects to detect this unique imprint of tidal dominance in making OOC comets 
discernable at the present epoch. 

Kaib and Quinn (2010) have found that w 5% of discernable OOC comets should have 60 < x, the 
majority of which have evolved from a primordial location in the IOC. It is unlikely that these comets have 
had their prior orbit perihelia outside the loss cylinder and therefore would not satisfy |AH| > |AH|^^j^. 
Taking account of expected errors in the binding energies we adopt 30 < a; < 60 as the defining interval for 
intermediate-^ comets. 

These clear signatures of galactic tidal dominance in making OOC comets discernable do not mean 
that we must abandon hope for detecting any impulsive imprint on the distributions, as we describe in Sec 
3. The characteristic S = ~1 continues to dominate for observed IOC comets, suggesting that they are 
likely to be daughters of OOC comets. S — —1 also remains dominant for other quality class OOC data 
(IB, 2A, 2B). However the dominance no longer correlates with the energy range 30 < x. We interpret this 
as evidence of significant errors in the determination of A for these quality classes, a conclusion made clear 
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by the excess of nominally unbound comets. 



2.3. Comparisons with recent modeling 

Kaib and Quinn (2009) have done the most complete modeling of the production of observable 
long-period comets (LPC) from the Oort cloud in an attempt to infer the Sun's birth environment. They 
closely follow whether an observed OOC comet was born in the OOC or in the IOC and find that roughly 
comparable numbers of each primordial population are made discernable after a suitable time lapse. 
Their binding energy distribution for the two discernable populations (their Fig. S4) differs in that the 
distribution for comets born in the IOC peaks at a; « 37 while the peak for comets born in the OOC peaks 
at a; « 27. They conclude that the IOC pathway provides "an important, if not dominant, source of known 
LPCs". The primordial origin of discernable comets is not important in our loss-cylinder modeling. The 
17th Catalogue indicates that only 2 of the 53 observed class lA comets with 100 < a: < 1000 have future 
orbits in the OOC, while 14 of 102 observed class lA OOC comets that have future orbits in the OOC. 
Therefore very few observed comets with 30 < a; < 60 are likely to have prior perihelia within the loss 
cylinder. We have used their Fig. S4 as a guide to the energy range most likely to show evidence of a weak 
impulsive component of the observed OOC. 

Rickman et al. (2008) have also embarked on an ambitious modeling program of the long term evolution 
of the Oort cloud which emphasizes a fundamental role for stellar perturbations. They demonstrate that 
over long timescales stellar impulses are needed to replenish the phase space of the OOC which is capable 
of being made discernable by the galactic tide. Massive star impulses serve to efficiently refill this feeding 
zone for a period of several 100 Myr and therefore provide one aspect of the synergy with the galactic tide 
that make these comets discernable. 

They also assert that "treating comet injections from the Oort cloud in the contemporary Solar System 
as a result of the Galactic tide alone is not a viable idea" . This statement follows from their observations 
that a tide-alone model evolving over the lifetime of the Solar System differs significantly at the present 
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epoch from a combined impulse-tide model. 

But their modeling does not shed light on the question of the dominant dynamical mechanism 
responsible for injecting OOC comets into the inner planetary region at the present epoch. Their work makes 
detailed predictions using results averaged over 170 Myr. Such a large time window inevitably includes 
many individual stellar perturbations which contribute directly to the time-averaged flux. Separating out 
the dominant perturbation making these comets discernable at the present epoch would be best obtained if 
they (i) "turn off' both perturbations at the present epoch, (ii) wait for the discerned flux to dissipate, and 
then (iii) alternately determine subsequent distributions produced by each perturbation separately turned 
on. This analysis was not performed. 



Indeed, one can compare our predictions (jMatese and Lissauer 



2004r )) for discernable distributions in 



A and the ^-dependent discernable zone refill efficiency with Rickman et al. (2008) (where they use the 
term "filling factor" in the same context as our "efficiency"). The tide-alone results at the beginning of 
their simulations (see their Fig. 7 and Table III) provide the most appropriate comparison since their phase 
space then is most nearly randomized and comparable to that used in our previous work. One finds that 
these two sets of distributions are in good agreement. 

Their combined modeling of the latitude distributions does predict a peak at | sini3| w 0.5, more nearly 
in agreement with observations shown in Fig. 2 (however they assert that their predictions do not agree 
with observations without any discussion of the nature of the discrepancies perceived) . If this difference 
with our tidal model prediction of a peak at | sini?| w 0.7 is supported in the future it may provide evidence 
that the in situ phase space refilling by stellar perturbations is indeed detectably incomplete at the present 
epoch, a point emphasized by them. This would not contradict our assertion that the unambiguous evidence 
is that the galactic tide dominates in making OOC come ts discernable at the present epoch . It remains 



2002 . 



20041 )') is no longer 



for them to demonstrate that the modeling adopted here (jMatese and Lissaueii 
viable in describing observations at the present epoch. We conclude our remarks by noting that neither of 
the above analyses have tracked the orbital characteristic S in the production of discernable OOC, a major 
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measure of galactic dominance in our work. 



3. Dynamics of a weak impulsive perturbation 
3.1. Theory 

The dynamics of a weak perturbation of a near-parabolic comet by a solar companion or field object is 
now considered. The change in the comet's specific heliocentric angular momentum induced by a perturber 
is given by 

1 1 



Hpert=A^p(Rxr)(^^^--j (3) 

where = GMp is the perturber mass located at heliocentric position r and R is the heliocentric comet 
position. In terms of the perturber true anomaly, /, we have r — + Cp cos /), with p — a(l — Cp^). We 
then obtain for near-parabolic comets 



where b = a^/l — ep^, the semiminor axis of the perturber. One then constructs 

AHpert = VA*0« — Q X I, (5) 
where the dimensionless integral I is taken around an interval of the perturber true anomaly A/ = /q — /; that 



corresponds to the comet orbital time interval between to = 0, the present epoch, and U = —1'k^A^I\iq^ 

The first term in I corresponds to the perturber interaction with the comet, the second with the Sun. For 
a specified perturber orbital ellipse, I is determined by the cometary Q and A as well as the companion's 
present value of the true anomaly, /q. Equation |3]is more convenient to use in an impulse approximation 
as has been done in a previo us analysis of the combined tide-stellar impulse interaction with the OOC 



(Matese and Lissauer 



20021 ) ) . Equation |4] is more appropriate if one wishes to include the slow "reflex" 



effects of a bound perturber on the Sun. 
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3.2. Combined tidal and impulsive interaction 

The galactic tidal perturbation and any putative point source perturbation of the Sun/OOC 
are, in nature, superposed in the course of a cometary orbit. For weak perturbations the two effects 
can be superposed in a vector sense. Hobs = Hprior + AHpert + AHtide = Hprior + AH„et- The 
cometary phase space of comets will have the prior loss cylinder distribution displaced by AHnet- 
The standard step function for the prior distribution of angular momentum is changed to a uniformly 
displaced distribution, similar to that illustrated in Fig. 1, but with AHtide replaced by AHnet, 

Qprior — Qobs = (^Hnet^ — 2Hobs ' ■^Hnet) f^fJ-Q- 



3.3. Weak impulsive effects on the discernable energy and spatial distributions 

Matese and Lissauer (2002) modeled the time-dependent changes in discernable OOC comet orbital 
element distributions that resulted from a weak stellar impulse. In particular, the in situ energy distribution 
was taken to be essentially unchanged, as described above. The number of comets in the largc-^ {x < 30) 
interval that became discernable after an impulse was found to be essentially unchanged from the case with 
no impulse, although the specific comets made discernable were changed. That is, the « 100% efficiency 
of refilling the discernable zone remains w 100% throughout the stellar impulse for large-A comets. This 
can be visualized in Fig. 1 when we consider a large AHtide that completely refills the loss cylinder and 
the comparative case where the large AHtide is modified by a weak AHpert- Independent of whether the 
impulse slightly increased or decreased the perturbation, a large | AHnet] will still tend to completely refill 
the loss cylinder, albeit with different comets. 

For the intermediate- A population the discussion is more subtle. Suppose that for a specific Q a comet 
population with a; « 35 will partially fill the discernable zone due to the tide (see Fig. 1). If some of those 
comets experience an impulse that increases | AHnet I, the discernable zone will be more completely filled. 
However other comets will be impulsed such that | AHnet] is decreased. This will in turn decrease the 
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number of discernable comets for this value of x leaving the efficiency for this x only moderately increased. 

Consider now the case for a population with x sa 45. The tidal perturbation will be smaller by a factor 
of ~ 0.4 from the x w 35 population, which may be inadequate to make any of these comets discernable 
in our loss cylinder model. In this case an impulsive torque preferentially opposed to the tidal torque will 
have no effect on the number of observed comets for this x since none would have been observed in its 
absence. But for those comets which have a weak impulsive torque preferentially aiding the tide, some will 
be made discernable that would not have been in the absence of an impulse. Therefore the efficiency for 
this X will increase from zero, and it will preferentially have S = —1. The net effect for a weak impulse is 
to create an enhanced observed OOC comet population along the track of the perturber that preferentially 
has intermediat e- A and 5' = — 1. Th e se feat ures have been demonstrated in detailed modeling for a weak 



stellar impulse (|Matese and Lissauer 



2002f) ). Rickman et al. (2008) also discuss in detail this aspect of 



synergy but do not consider the importance of the characteristic, 5, in this discussion. 

In reality, the step function distribution of prior orbits for OOC comets used in the loss cylinder 
model is a crude first approximation. A small fraction of original OOC comet orbits recede from the 
planetary region as future OOC comets with perihelia inside the loss cylinder. This does not obviate the 
arguments invoked above. Along with the original energy errors associated with observational uncertainties 
and outgassing effects, we can understand why the observed spread in original energies seen in Fig. 3 is 
somewhat larger t han pr edicted in our loss cylinder model and is consistent with more realistic modeling 



(jKaib and Ouinnl (| 2009I )). 



4. Observational evidence for an impulse 

Matese et al. (1999) noted an excess of OOC comet major axes along a great circle roughly centered on 
the galactic longitudinal bins 135° and 315° using 11th Catalogue data of quality classes lA+lB. Further, 
they found that the excess was predominantly in the intermediate-A population and had a larger proportion 
of 5* = — 1, consistent with it being impulsively produced. In Fig. 4 we now display the distribution of the 
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aphelia longitudes for our 17th Catalogue data of quality class lA. We see that the excess remains in the 
present data. 



4.1. Inferring the orbit orientation of a putative bound perturber 

Our first investigation relates to the orientation of the overpopulated band. There is no obvious a 

priori reason why an impulsive track should pass through the galactic poles. Further, choosing longitude 
bins for comparison has an obvious bias in that it preferentially excludes solid angles near the poles where 
the tide is weak. 

We therefore perform a great circle fit that counts the number of major axes within an annular band of 
width ±9.6°, which has the same solid angle, 4n/6, as the two overpopulated aphelia longitude bins in Fig. 
4. The data that we now analyze includes the 35 (of 102) 17th Catalogue comets listed in the Appendix 
which are intermediate- ^4 comets, 30 < a; < 60, with S = —1. As described above, this is the subset most 
likely to exhibit evidence of a weak impulsive component. Of these 35 comets, 15 were included in the 11th 
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Fig. 4. — The distribution of aphelia longitude, L, for class lA OOC comets. 
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Catalogue considered by Matese et al (1999). Counting is done for all possible great circle orientations 
with normal directions stepped in a grid covering solid angles (5°^ • The process is repeated for the 11th 
Catalogue data, the complete 17th Catalogue data as well as for the 20 comets in the 17th — 11th Cataogue 
subset. Great circle orientations are denoted by the galactic longitude of ascending node, Cl, and the 
galactic inclination, i. A similar plot holds for opposite great circle normal vectors 180° — i with O — 180°. 
11th catalogue, n>10 17th-llth catalogue, n>10 17th catalogue, n>20 




300 310 320 330 340 300 310 320 330 340 300 310 320 330 340 

n n n 



Fig. 5. — Orientation angles (O =galactic longitude of ascending node, i= galactic inclination) of great circle 
bands of width ±9.6° which include 10(20) or more major axes of intermediate-A comets with S = —1 . 

For the 11th Catalogue data we can find a great circle with as many as 13 (of 15) axes inside the band. 
We show all orientations which include 10-13 major axes in Fig. 5. For the 17th — 11th data we show all 
orientations that include 10-11 (of 20) major axes. One observes that the peak overpopulated great circle 
normal direction occurs in the same general area of the hemisphere for both the 11th and 17th — 11th 
data. The probability that these solid angles would overlap if the two sets of data were uncorrelated can be 
estimated to be w 0.02. For the 17th Catalogue data we show those directions that include 20-22 (of 35) 
major axes. The statistical significance of these observations is discussed below. 

The best fit great circle normal direction with a ±9.6° band containing 22 axes of the 17th Catalogue 
data is centered at z = 103°, fl = 319° (or its opposite direction). A measure of the uncertainty in the 
great circle orientation can be inferred from Fig. 5. In ecliptic coordinates the best fit great circle plane is 
specified by = 133°, fie = 190°, or its opposite direction. 
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In the Appendix we list the 17th Catalogue data for the 102 class lA comets that constitute our 
complete set. For each comet we include the binding energy parameter, x, the tidal characteristic, S, the 
perihelion distance, q, and the galactic angular orbital elements B, L and a. The last column gives the 
magnitude of the angular separation of the cometary major axis from the best fit great circle plane, 7. 
Comets having angular separations within ±7 fall inside a band of solid angle 47r| sin7|. Comets with year 
designation prior to 1995 constitute our 11th Catalogue data. 



Fig. 6. — The scatter distribution for aphelia directions of 17th Catalogue class lA OOC comets having a 
binding energy parameter, x = 10^ AU/A, in the interval 30 < x < 60 and tidal characteristic S=-l. Solid 
curve: ecliptic plane. Dashed curve: best fit perturber plane with i = 103°, O = 319°. Black dots: aphelia 
directions within a band of width ±9.6° covering 1/6 of the celestial sphere and centered on the best fit great 
circle path. Gray dots: exterior to the band. 

In Fig. 6 we illustrate the scatter of aphelia directions in galactic coordinates with the best fit great 
circle arc {i = 103°, O = 319°) shown. The data illustrated is for the 17th Catalogue with 30 < a; < 60 and 
S = —1. Major axes are distinguished if they are within the ±9.6° band. The distribution of major axes 
around the great circle arc indicates comparable numbers in each of the 4 quadrants. This is inconsistent 
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with an overpopu lation produced by a weak s tellar perturbation which produces an overpopulated arc 



< 150° in length (jMatese and Lissauer 



20021 ) V In Fig. 7 we show the cumulative number of major axes 
interior to the solid angle, 47rsin7, subtended by a band of width ±7 centered on the best fit great circle. 
For comparison, results are shown for the 17th Catalogue data as well as the corresponding 11th Catalogue 
data. 

Using an approach that ignores the characteristic S, Murray (1999) has argued for a solar companion 
that is the source of a impulsed observed OOC population with 40 < x < 100 whose axes are concentrated 
along a great circle that is different from that found here. He found that all class lA comets in the IQth 
Cometary Catalogue had axes within a ±30° band. His point is essentially that all of these observed comets 
must have been injected into the loss cylinder by an impulse of a bound companion. We have taken the 17th 




Sin 7 



Fig. 7. — The cumulative distribution of the number of major axes of class lA OOC comets having 
30 < a; < 60 and S — ~1 that fall within an annular band of solid angle 47rsin7 centered on the best 
fit great circle i = 103°, = 319°. Black: 17th Catalogue data. Gray: 11th Catalogue data. Solid lines 
are for random distributions. Dashed line indicates annular band of width ±9.6° which includes 1/6 of the 
celestial sphere. 
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Catalogue data and considered all lA (40 < x < 100) comets that also were included in Murray's paper 
which analyzed IQth Catalogue data. A fit to the original data is found to agree with Murray's result, a best 
fit companion orbit plane with a galactic inclination oi i — 37° and galactic ascending node of = 175° 
that includes all 12 IQth Catalogue axes within a band of ±30° which encloses 1/2 of the celestial sphere. 
When we repeat the fit with the same band orientation and width we find that 10 of 25 additional comets 
in the llth Catalogue are enclosed. The concentration originally perceived is seen to have disappeared. 



4.2. IR observational limits on the companion's mass and distance 

The absence of detection in the IRAS and 2mass data bases place limits on the range of possible 
masses and present distances of the proposed companion. For relevant values the 2mass limits are not very 
restrictive. Ass uming a J-band detect ion magnitude limit of 16 and theoretical models of a 4.5 Gyr old 



2-10 Mj object ([Burrows et al. 



20031) ). we find that the mass would have to be greater than 7 (lO)Mj and 



closer than 6000 (25, 000) AU for a possible 2mass identification. 

For Jovian masses, the strongest current limits on distance are found from the IRAS Point Source 
and Faint Source Catalogs. We have performed a non-parallax search of both IRAS catalogs using the 



VizieR Service (http://vizier.u-strasbg.fr/cgi-bin/VizieRl. The search was based on the theoretical models 



of Burrows et al. (2003). Search criteria included rejection of previously identified point sources, a PSC 12 
/im flux greater than 1 Jy and a FSC fiux greater than 0.3 Jy. The latter value is based on a completion 
limit of 100% (Moshir et al. 1992). Only flux quality class 1 sources were included. We further required 
that there be no bright optical or 2mass candidates within the 4cr error limits and that the absolute value 
of the galactic latitude be greater than 10 degrees, as required for sources to appear in the FSC. Distant 
IR galaxies were excluded by requiring the 25 /im flux to be greater than the 60 /im fiux. The PSC search 
yielded a single candidate, 07144-1-5206, with galactic coordinates B — 25° and L — 165°, not far from the 
great circle band. The 12 and 25 fim band fiux (0.8 Jy and 0.5 Jy) of this source is consistent with a « 6 
Mj mass at « 6000 AU. However, the FSC, which is the definitive catalog for faint sources, associated this 
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PSC source with another source located 80 arc sec away. This position corresponded to an extremely bright 
(J=6) 2mass source. The negative IRAS search results suggest that an object of mass 2 (5) Mj must have 
a current distance r > 2000 (10,000) AU, respectively. 

4.3. Dynamical inferences of other perturber properties 

Although the orbit normal of the putative companion is tightly constrained, other properties are less 
so. If it exists, a solar companion most likely formed in a wide-binary orbit. The near-uniform distribution 
of the overpopulation around the great circle suggests that any putative companion is likely to have a 
present orbit that is more nearly circular than parabolic {Cp^ < 0.5). The near-circular implication cannot 
remain true over the solar system lifetime since the eccentricity and inclination osculate significantly due to 
the tide. Further, the semimajor axis will be affected by stellar impulses over these timescales. 

The implication that Cp^ < 0.5, is in fact, consistent with the present galactic inclination of the 
perturber inferred here, i w 103°. This follows from the near-conservation of the z component of galactic 
angular momentum of objects in the IOC and OOC in the intervals between strong stellar impulses, 
Hz oc b cosi « constant. Since the present value of | cos 103° | « 0.22 is at the low end of a random 
distribution, < | cosi| < 1, a larger primordial value of | cosi| implies that the present value of Cp is 
reduced from its primordial value. 

In Fig. 8 we show a representative time dependent companion orbit assuming present {to = 0) galactic 
orbital elements i = 103°, fl = 319°, ui = 0, a = 15,000AU and Cp = 0.5. We integrate back in time 
assuming an unchanging galactic environment and neglecting impulsive perturbations, both of which cannot 
be ignored for more reliable results. The purpose of the calculation is to demonstrate that significant 
changes in the eccentricity and inclination occur over the solar system lifetime for a companion with 
a > 10^ AU. 

There is no evidence of overpopulation within the perceived band of OOC or IOC comets with 60 < x 
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where |AHtide| < lAHminl- Here we must distinguish two subpopulations of discernable comets with 
60 < a; (assuming negligible errors in cc), (i) comets that have been directly injected from beyond the loss 
cylinder by the perturbation and are making their first orbit as discernable comets, or, (ii) comets that are 





Fig. 8. — Integrating back in time for 4.6 Gy the osculating orbit of a putative companion with present {t = 0) 
galactic orbital elements i{0) = 103°, 17(0) = 319°, w(0) = 0, a = 15,000AU, ep(0) = 0.5, g(0) = 7500 AU 
assuming unchanging galactic environment and neglecting impulsive perturbations. Shown are the perihelion 
distance, q{t), and the semiminor axis, b{t), along with the galactic and ecliptic inclinations, i{t) and ie{t). 
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daughters of a; < 60 comets which had a previous perihehon passage well inside the loss cylinder. The latter 
is the commonly accept ed pathway for producing disce rnable IOC comets, but requires an ad hoc imposition 



19991 )) to explain the relative numbers of observed IOC 



of a fading mechanism (jWiegert and Tremaine 
and OOC comets. 

For a comet major axis inclined to the perturber orbit plane by an angle 7 |AHpert|max varies 
approximately as I/7, and the cross section for scattering comets varies approximately as 7^. Thus one 
expects that a perturber in a crossing orbit with some 60 < a; comets would produce a narrower band of 
detectable flux perturbed with |AHpert| > |AHmin|- The absence of such a directly impulsed population 
(i) is not inconsistent with the notion that the depletion rate of this population in crossing orbits within the 
great circle band by the perturber would be larger than the rate of refilling of the band phase space by the 
galactic tide. 

More difficult to explain is the absence of an enhanced daughter population (ii) in the perceived 
band. A possible explanation is that the reduced impulse/reflex interaction can remove daughters by 
a random-walk of their perihelia from inside the discernable zone leading to ejection over one or more 
subsequent orbits. A direct numerical integration, including a putative companion, of the type done by 
Kaib and Quinn (2009) is needed to clarify this matter. This scenario would require a to be at the small 
end of the allowed range. 

Setting the critical perturbation for 7 = 9.6° equal to the minimum change required to enter the 
discernable zone, Eq. [5] yields 

lAHpertlcrit = ^/Jjk^a — |Q X I|crit = |AHmi„| ~ Vmo 5.4AU. (7) 

This provides a means of estimating the perturber properties needed to weakly impulse the OOC comet 
population and assist the tide in making a comet with 30 < a; < 60 discernable. Here |Q x Ijcrit is taken to 
be ~ the peak values of |Q x I| sampled over all orientaions of Q that are inclined to the great circle by an 
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angle 7 = ±9.6° 
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Fig. 9. — An example case showing the single-orbit impulse strength |Q x I| for comets with A = 22,000 
AU and Q inclined to the perturber plane by 7 = 9.6°. The positions a = 0(180°) correspond to comet 
orientation Q closest to the perturber perihelion (aphelion). Perturber parameters are a = 15,000 AU, 
Cp = 0.5, and /o = 215°, the present true anomaly of the perturber. 

In Fig. 9 we show an illustrative plot of |Q x I| for perturber parameters a = 15,000AU, Cp = 
0.5, fo = 215° and comet A = 22, OOOAU. A set of major axis orientations, Q, are considered which are 
all inclined along a conical surface at an angle 7 = 9.6° to the perturber plane. The specific orientation 
(J = 180° corresponds to the case where Q is closest to the perturber aphelion direction. The arbitrarily 
chosen present location of the companion is fo = 215° so that it has recently passed aphelion. Two peaks 
are shown corresponding to the two comet axes Q that are most strongly perturbed on the way out and 
on the way in during the prior cometary orbit. The peak at u = 195° locates the perturber position when 
comets were perturbed on their inward leg. The ratio of the comet/perturber periods is f»1.8 in this case 
and the second peak at cr = 35°, corresponds to the outward leg when the perturber was on its previous 
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orbit. The importance of the reflex term can be gauged from the background contribution. In general the 
peak corresponding to the inward comet leg will slightly lag the present perturber position if the perturber 
is presently more nearly at aphelion. The second peak corresponding to the outward leg will, in general, 
tend to be randomly positioned relative to the inward peak ii A > a. Each peak contains comet axes Q 
that have angular momentum component impulses that both aid the tide and oppose the tide. 

The results shown are typical for ah Cp < 0.5, 90° < fo < 270° when 10,000AU <a< 30,000AU. The 
exceptions include the less likely cases when the perturbations on both the outward and inward legs of a 
comet orbit find the perturber at the same general location / ~ (t creating a larger net impulse. We shall 
adopt a typical value of |Q x Ijcrit ~ 6 — 12 for all possible orbit parameters above. Smaller values in this 
range typically occur for smaller values of a which have higher relative velocities at orbit crossing. Another 
aspect of the variation is that the small refiex term can add or subtract from the direct impulse term. 

Inserting Mp = 5Mj into Eq. |S] we obtain a < 6000 AU. Assuming e^^ < 0.5, this is marginally 
inconsistent with IRAS observational limits. Smaller masses should not confiict with IRAS. Thus we 
adopt an approximate range of perturber parameters IMj < Mp < 2Mj for a « 30,000AU extending to 
2Mj < Mp < 4Mj for a « 10, OOOAU. This compaii i on ma s s range is consistent w i th the minimum Jeans 



"1 r 



mass (1-7) Mj as variously calculated (|Whitworthl (| 20061 ). 



Low and Lvnden-Belll (| 19761) ) 



A solar companion remains a viable option (iMatese et al. 



19991) 



Horner and Evans 



20021 )). In 



addition, we find that an 4 Mj companion in an orbit such as that shown in Fig. 8 would be capable 
of adiabat ically detaching a sca t tered disk EKBO, pr oducing an object with orbital characteristics similar 



to Sedna (Matese et al. 



20061 ): 



Gomes et al. 



20061 )). Allowing for the possibility that the perturber 
was more tightly bound primordially, smaller masses are then allo wed. It may a l so be possible to explain 



the misalignment of the invariable plane with the solar spin axis ([Gomes et 



20061 )) if the putative 



companion was more t ightly bound primordially. Such a companion would have a tempera ture of ~ 200 K 



(Burrows et al. 



20031 )) and will be easily seen by the Wide-field Infrared Survey Explorer ([Wright 



20071) ) 



recently launched. Gaia, which will use astrometric microlensing, may also be able to detect the putative 
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companion and will be launched in 2011 (|Gaudi and BloomI () 2005D ) 



Matese and Lissauer (2002) investigated the frequency of weak stellar showers as well as the patterns 
of discernable OOC comet orbital elements. They found that stellar showers that produce a > 20% peak 
enhancement in the background tidal flux occur with a frequency of approximately once every 15 Myr. The 
half-maximum duration is ~ 2 Myr. Thus it is only moderately unlikely that we are presently in a weak 
stellar shower of this magnitude. 

Flux enhancements of this magnitude were found to extend over an angular arc ~ 150° and have a 
full-width half-max of « 50°. Weak stellar impulses never extend over a larger arc. More simply put, to get 
the same stellar-induced flux enhancement as inferred from the data here, the enhanced region will have 
an extent that is « half as long and « twice as wide as that observed. Of course one must consider the 
possible alignment of two weaker coincident stellar impulses (or a single weaker statistical anomaly and a 
coincident single weaker stellar impulse) that happen to line up on opposite hemispheres. These will be 
improbable. For example, suppose we assume that there will always be two weak stellar impulses enhancing 
the observed background tidal flux, each with half the observed enhancement. The probability that both 
stellar orbit planes will be aligned to within ±9.6° is 0.014. 



4.4. A Monte Carlo test for the likelihood of the conjecture. 



For e ach of the three data sets d iscussed above we perform a variation of the Monte Carlo test done 
elsewhere 



Horner and Evans 



20021) ) which considered the original conjecture. Each comet in the data 
maintains its measured value of B, but is randomly assigned a value for L. The best fit analysis is performed 
to find the orientation of the great circle band of width ±9.6° which maximizes the number of major 
axes interior to the band, just as was done for the real data. We then record that number, n. Repeating 
the process for each data set 5000 x we can find the Monte Carlo probability that > n major axes would 
fall in the band if, in fact, L was truly random. Results are given in Table 1. As an example, for the 
17th Catalogue data, when we randomize L for the 35 comets we found an unconstrained (i.e. randomly 
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oriented) band of width ±9.6° which maximized the number of major axes contained in the band at 18 or 
more in 30 cases of the 5000 tested. 

Table 1: Monte Carlo probability ^Pmc{^ n) for Catalogue data set £ 



£ 
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11 


12 


13 


nth 


15 


0.018 


0.0020 


0.0002 


<0.0002 


nth -nth 


20 


0.24 


0.078 


0.015 


0.0014 




n 


18 


19 


20 


21 


mil 


?,r> 


O.OOGO 


0.0018 


0.000<S 


<0.0002 



A preliminary discussion of the persistence of the enhancement is as follows. The Monte Carlo 
probability, ^^*^'Pmc{'> 10) = 0.018 for the original data set is superficially significant. The solid angle 
over which this occurs can be seen in Fig. 5 to be < 2% of the celestial sphere. As many as n = 13 of 15 
axes can be found in a band for the actual data. Although the Monte Carlo analysis for the nth — nth 
Catalogue data set is, in itself, not significant since at most 11 comets appear in a band for the actual data 
and i^*''-i"''Pmc(> 10) = 0.24, the fact that the cone of normal directions that maximizes the counts for 
this data set overlaps with the cone that maximizes the nth Catalogue data set is significant. To properly 
gauge this significance we discuss a Bayesian analysis of the likelihood of the present conjecture. 

4.5. Bayesian inference 

The anomalous comet data discussed above can be explained by the putative companion. We now 
wish to compare this hypothesis with other hypotheses that might also explain the data. The most obvious 
alternative explanation is that the anomaly is entirely a statistical fluke. The other dynamical explanation 
is that of a weak stellar shower assisted by the galactic tide. Arguments against this possibility have been 
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given earlier and in the following analysis we assume the probability of a shower explanation is negligible. 



We consider the hypothesis, Ji, that the there exists a Jovian mass solar companion orbiting in the 
Oort cloud with parameters in the previously discussed ranges that are allowed on the basis of IRAS/2mass 
observations. The null hypothesis, -^H, is that no such object exists. £ constitutes the evidence considered, 
the 30 < a; < 60, S = —1 data sets described in Fig. 5 and Table 1. The a priori relative probability of 
hypothesis H, before we consider evidence £, is symbolically defined as V{H) /'P{^H). 

IRAS and 2mass observational constraints prohibit present locations r < 2000 : 10,000 : 25,000 AU 
for Mp = 2 : 5 : 10 Mj respectively. Dynamical constraints previously discussed suggest a companion 
mass in the range Mp sa 1 - 4 Mj and semimajor axis in the range a « 30, 000 — 10, 000 AU, respectively. 
Qualitatively 'P{H)/'P{-^H) is the ratio of the specified companion parameter space area {Mp,a) to the 
complementary, non-IRAS(2mass)-excluded parameter space area. Next we try to roughly estimate this 
from the non-Catalogue evidence. 

We assume that there is a certainty that the sun has an as-yet-to-be-discovered maximum mass 
companion with mass between Mpiuto ■^ — ^ 20Mj, the upper limit being fixed by IRAS/2mass observations. 
Zuckerman and Song (2008) give a mass distribution for brown dwarfs down to 13 Mj of dN/dM cx M^^-"^ . 
It m ight be assumed th at this power law roughly holds down to the minimum Jeans mass limit of 1 - 4 



Mj (jWhitworth 



20061 )) but the mass distribution below 1 Mj is completely unknown. For simplicity we 



assume dN j dM oa M ^ over the entire mass range. 

Next we estimate the probability distribution in semimajor axis for this companion, assumed to be 
in the interval 10^ < — > 10^ AU. We base this estimate on the known distribution of semimajor axes of 
wide binary stars. At least 2/3 of solar type stars in the field reside in binary or multiple star systems 

)). The fraction of these s y stems which have semimajor a xes a > 1000 



(jDuquennov and Mavor 



199 



AU (wide binaries) is « 15% (|Kouwenhoven et 



20101) 



Duquennov and Mavor 



19911) ), and includes 



periods up to 30 Myr. The origin of these wide binaries may be capture during cluster dissolution 



(Kouwenhoven et al 



20101 )). For large semimajor axes the observed distribution is well approximated by 
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{dN/ da) oc a ^ and we adopt this power law for all masses over the entire interval of a. 

The simplifying assumptions we have made are equivalent to assuming a uniform probability 
distribution over the entire parameter space (Mp, a) under consideration 

— — — = constant. 

d log Mp d log a 

Thus the combined probability that the sun has a companion of mass 1-4 Mj with semimajor axis 
30,000 - 10,000 AU is simply the ratio of this parameter area to the entire area after the IRAS constraint 
has been utilized. In the following Bayesian analysis we use this value for the a priori probability ratio, 

p{n)/p{-^n) « 0.01. 



The likelihood ratio of Ti given £ is 



where V{£\'H) is the conditional probability of evidence £ if hypothesis V. is true. More meaningful is the a 
posteriori relative probability, a product of the a priori relative probability and the likelihood ratio 

^^^1') -Ah.:^. (10) 



V{-^H\£) "-^ v{-^Hy 

4-5.1. Analysis of the 17th Catalogue alone 
In this case we adopt the Monte Carlo results of Table 1 for the likelihood ratio 

A = 1 

n:17th i7t/.p^^(> n) 

to obtain an a posteriori probability ratio 

v(n\mh) 0.01 



Inserting the results from Table 1 the a posteriori relative probability ratio is « 10 for n = 20 with 
corresponding changes for other values of n. 
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4-5.2. Separate analyses of the 11th and 17th— 11th subsets 

For £ = 11th we make the same assumption about the a priori relative probabihty to obtain 

V{H\llth) _ 0.01 
V{--H\llth) ~ i"'*Pmc(> n) 

from which we obtain « 0.5 for n = 10 with corresponding changes for other values of n. 

For £ = nth- 11th we must take into account that the 11th Catalogue data defined a specific cone of 
directions in space. In this case we adjust our likelihood ratio estimate of the Monte Carlo results of Table 
1 by the probability that the cones of the orbit normals containing maximum counts would overlap which 
is Ri 0.02 to obtain an a posteriori probability ratio 

V{H\mh-llth) _ 0.01 
r{-nn\17th-Uth) ~ 0.02 [ l7t'^-llf^pJ^^c,(> „)] 

from which we obtain results of w 2 for n = 10 with corresponding changes for other values of n. 

An alternative approach is to allow our 11th Catalogue inference to adjust the a priori relative 
probability for £ = 17th — 11th {if the cones overlap) to equal the a posteriori result for £ = 11th 

V{H\mh-llth) 05 

V{--H\17th-llth) ~ i7t/j-iif?»p^^(> n) 

which also yields w 2 for n = 10. 

All of these inferences can be scaled if a value of the a priori relative probability is chosen different 
from 0.01 based on other evidence than considered here (IRAS/2mass observational limits and stellar/brown 
dwarf mass and semimajor axis distributions). Specifically this value might be increased if one believed 
that the Jovian companion hypothesis, "H, provides a plausible explanation for producing Sedna, and it 
would be substantially reduced if one believed that the hypothesis is incompatible with the absence of an 
enhanced IOC daughter population along the perceived arc. 

We can estimate the size of the impulsive enhancement in the present data for 102 class lA comets as 
follows. Although only 13 of 35 comets lie outside of the ±9.6° band for our best fit, statistically significant 
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results obtain if as many as 17 of 35 were outside the band. So, conservatively, we would expect w 3 inside 
the band if there were no impulsively enhanced component. The excess of « 15 — 19 inside the band is 
« 20% of the complementary number of 87-83 for the complete lA data listed in the Appendix. The same 
number can be obtained by extrapolating back to zero the large sin 7 portion of the cumulative curve in 
Fig. 6. 

5. Summeiry and conclusions 

We have described how the dynamics of a dominant galactic tidal interaction, weakly aided by an 
impulsive perturbation, predicts specific properties for observed distributions of the galactic orbital elements 
of outer Oort cloud comets. These subtle predictions have been found to be manifest in high-quality 
observational data at statistically significant levels, suggesting that the observed OOC comet population 
contains an « 20% impulsively produced excess. The extent of the enhanced arc is inconsistent with a weak 
stellar impulse, but is consistent with a Jovian mass solar companion orbiting in the OOC. A putative 
companion with these properties may also be capable of producing detached Kuiper Belt objects such as 
Sedna and has been given the name Tyche. Tyche could have significantly depleted the inner Oort cloud 
over the solar system lifetime requiring a corresponding increase in the inferred primordial Oort cloud 
population. A substantive difficulty with the Tyche conjecture is the absence of a corresponding excess in 
the presumed IOC daughter population. 
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Appendix 



Table 2: Class lA OOC orbital elements. All angles are in degrees. 



Comet name 


X 


S 


g(AU) 


B 




L 




a 


7 


C/1996 Jl-A(B) 


-510(-1) 


-1 


1, 


.2978 


15, 


.9 


202 


.6 


55, 


.0 


51.13 


C/1996 El 


-42 


1 


1, 


.3590 


-22, 


.6 


304, 


.3 


319, 


,2 


18.39 


C/1942 C2 


-34 


1 


4, 


.1134 


-24, 


.5 


18, 


.8 


335, 


.2 


42.27 


C/1946 CI 


-13 


1 


1, 


.7241 


-60, 


.3 


313 


.1 


276, 


.7 


14.17 


C/2005 Bl 


-13 


1 


3, 


.2049 


-17, 


.6 


82, 


.6 


12, 


,7 


38.28 


C/1946 Ul 


-1 


1 


2, 


.4077 


49, 


.0 


347 


.2 


255, 


.0 


28.13 


C/1997 BA6 





-1 


3 


.4364 


27, 


.7 


123 


.0 


52, 


.9 


20.05 


C/1993 Ql 


3 


1 


0, 


.9673 


41, 


,2 


112 


.7 


174, 


,6 


28.22 


C/2005 Kl 


7 


1 


3 


.6929 


-17, 


.3 


227 


.4 


312, 


.1 


85.45 


C/2003 T3 


10 


1 


1, 


.4811 


25, 


.5 


293 


.6 


148, 


.9 


16.25 


C/2006 L2 


10 


-1 


1, 


.9939 


-50, 


.0 


225 


.1 


98, 


,2 


52.89 


C/1974 VI 


11 


-1 


6, 


.0189 


28. 


.0 


308 


.9 


63, 


.1 


2.61 


C/1988 Bl 


13 


-1 


5. 


.0308 


-47, 


.6 


316 


.3 


122, 


.1 


11.37 


C/2006 Kl 


13 


1 


4, 


.4255 


66, 


,1 


56, 


.2 


247, 


,3 


36.69 


C/2004 X3 


14 


1 


4, 


.4023 


-36, 


.8 


80, 


.9 


24, 


.1 


31.90 


C/2005 E2 


14 


1 


1, 


.5196 


42, 


.0 


312 


.6 


134, 


.6 


3.99 


C/2006 0F2 


15 


1 


2, 


.4314 


10, 


,5 


322 


.5 


133, 


,0 


5.75 
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Comet name 


X 


S 


?(AU) 


B 


L 


a 


7 


C/1991 F2 


16 


-1 


1.5177 


29.7 


90.0 


345.1 


48.65 


C/2005 Gl 


16 


-1 


4.9607 


-55.3 


297.9 


100.6 


22.60 


C/1916 Gl 


17 


-1 


1.6864 


17.9 


225.9 


30.1 


58.95 


C/1956 Fl-A 


17 


1 


4.4473 


-25.4 


180.0 


62.1 


42.33 


C/1944 K2 


18 


-1 


2.2259 


-31.8 


124.3 


159.3 


5.25 


C/2007 01 


18 


-1 


2.8767 


2.1 


197.3 


7.6 


55.06 


C/1935 Ql 


19 


1 


4.0434 


11.8 


248.0 


115.3 


58.81 


C/1999 J2 


19 


-1 


7.1098 


-49.9 


234.4 


259.3 


52.80 


C/1973 El 


20 


-1 


0.1424 


17.0 


207.3 


46.3 


53.15 


C/1984 W2 


20 


1 


4.0002 


33.1 


79.5 


101.0 


55.69 


C/2001 Gl 


20 


1 


5.1046 


-8.3 


135.0 


344.1 


2.02 


C/2004 PI 


20 


-1 


6.0141 


20.2 


208.7 


63.7 


51.26 


C/1922 Ul 


21 


-1 


2.2588 


16.3 


270.6 


52.4 


39.46 


C/1997 Al 


21 


1 


3.1572 


-19.9 


352.3 


306.7 


25.22 


C/2003 K4 


23 


1 


1.0236 


-46.8 


104.3 


307.2 


12.49 


C/2005 Ql 


23 


-1 


6.4084 


-8.8 


321.7 


258.1 


0.65 


C/1947 SI 


24 


1 


0.7481 


10.9 


185.4 


170.2 


40.61 


C/1978 HI 


24 


1 


1.1365 


-20.7 


124.4 


30.2 


8.65 


C/1999 K5 


24 


-1 


3.2554 


34.1 


129.8 


271.8 


14.81 


C/20G1 Gl 


24 


-1 


8.2356 


-46.8 


113.6 


118.9 


7.02 


C/2006 HW51 


25 


1 


2.2656 


-35.4 


162.3 


83.4 


26.39 


C/1914 Ml 


27 


-1 


3.7468 


-1.4 


200.9 


100.9 


59.81 


C/1980 El 


27 


1 


3.3639 


-18.2 


177.0 


29.9 


39.78 


C/2004 YJ35 


27 


-1 


1.7812 


-38.2 


346.3 


122.2 


12.23 



38 



Comet name 


X 


S 


g(AU) 


B 


L 


a 


7 


C/1992 Jl 


28 


1 


3.0070 


-43.8 


316.9 


0.3 


10.47 


C/1913 Yl 


29 


-1 


1.1045 


-52.4 


283.0 


237.4 


31.83 


C/1987 W3 


29 


-1 


3.3328 


64.7 


333.5 


62.5 


17.92 


C/2001 K5 


29 


1 


5.1843 


-30.3 


223.1 


78.8 


71.80 


C/1989 XI 


32 


-1 


0.3498 


-41.0 


325.6 


218.0 


3.57 


C/1978 Al 


33 


-1 


5.6064 


-31.4 


142.3 


167.8 


9.48 


C/1993 Kl 


33 


-1 


4.8493 


2.3 


130.9 


347.4 


8.38 


C/1948 El 


34 


1 


2.1071 


-38.6 


250.2 


83.2 


58.25 


C/1948 Tl 


34 


-1 


3.2611 


17.9 


324.0 


306.9 


8.63 


C/2006 Ql 


34 


-1 


2.7636 


-45.5 


112.2 


134.0 


8.48 


C/1925 Fl 


35 


1 


4.1808 


-40.9 


65.6 


327.8 


33.96 


C/2002 J4 


35 


-1 


3.6338 


32.4 


162.6 


342.6 


12.07 


C/1974 Fl 


36 


-1 


3.0115 


21.7 


140.3 


10.7 


3.59 


C/2003 S3 


36 


-1 


8.1294 


11.6 


345.3 


351.9 


27.95 


C/2004 Bl 


36 


1 


1.6019 


27.7 


167.1 


138.4 


17.55 


C/1950 Kl 


37 


-1 


2.5723 


-43.8 


137.7 


269.9 


8.02 


C/1999 U4 


37 


-1 


4.9153 


-30.2 


322.6 


183.8 


3.42 


C/2006 El 


37 


-1 


6.0406 


31.4 


140.5 


290.6 


21.24 


C/2006 PI 


37 


1 


0.1707 


-5.0 


342.2 


72.9 


5.47 


C/1999 Fl 


38 


-1 


5.7869 


15.2 


99.7 


341.0 


68.86 


C/1999 Ul 


38 


-1 


4.1376 


25.1 


67.4 


291.1 


40.91 


C/2000 CT54 


38 


-1 


3.1561 


36.7 


145.1 


8.0 


2.96 


C/2002 L9 


38 


-1 


7.0330 


51.6 


147.6 


292.5 


4.93 


C/1954 Yl 


39 


-1 


4.0769 


-23.5 


314.1 


232.3 


9.55 
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Comet name x S ^(AU) B L a 



C/1925 Gl 


40 


1 1.1095 


10.5 


238.3 


151.2 


64.77 


C/1997 J2 


40 


1 3.0511 


-0.5 


261.1 


332.6 


55.83 


C/1954 02 


42 


1 3.8699 


-18.4 


336.4 


53.3 


11.82 


C/1979 M3 


42 


1 4.6869 


14.1 


207.0 


126.6 


55.24 


C/2000 Kl 


42 


-1 6.2761 


-21.6 


190.7 


144.5 


52.50 


C/2006 S3 


42 


-1 5.1309 


-15.3 


188.0 


219.4 


50.21 


C/1946 PI 


44 


-1 1.1361 


-20.0 


126.6 


108.7 


6.85 


C/1999 Yl 


44 


-1 3.0912 


62.7 


289.2 


326.2 


1.31 


C/2000 Al 


44 


-1 9.7431 


-33.0 


19.7 


171.9 


36.22 


C/2005 EL173 


44 


-1 3.8863 


23.3 


50.4 


350.5 


79.60 


C/1932 M2 


45 


-1 2.3136 


-14.6 


146.5 


172.1 


10.33 


C/1987 HI 


46 


-1 5.4575 


-17.6 


190.9 


161.4 


53.02 


C/1888 Rl 


48 


-1 1.8149 


59.5 


313.7 


289.0 


8.54 


C/1983 01 


48 


-1 3.3179 


55.2 


324.4 


355.4 


13.71 


C/1973 Al 


49 


-1 2.5111 


-50.7 


106.2 


220.9 


9.24 


C/1989 Yl 


49 


1 1.5692 


-6.1 


328.9 


96.2 


8.22 


C/2007 Dl 


50 


-1 8.7936 


-41.3 


64.4 


157.8 


33.83 


C/2003 WT42 


51 


-1 5.1909 


-56.8 


353.7 


171.0 


6.59 


C/2000 01 


52 


-1 5.9217 


23.5 


313.8 


286.2 


25.23 


C/2000 SV74 


52 


-1 3.5417 


-21.5 


296.7 


193.5 


0.53 


C/1999 F2 


54 


-1 4.7188 


-46.0 


135.3 


112.9 


6.83 


C/1999 S2 


56 


-1 6.4662 


22.2 


153.2 


322.9 


7.83 


C/2006 K3 


56 


1 2.5015 


46.6 


22.7 


261.9 


49.79 



40 



Comet name x S g'(AU) B L a ^ 



C/1976 D2 


59 


1 6.8807 


-44.5 


120.9 


315.8 


3.36 


C/1987 Fl 


59 - 


1 3.6246 


-26.5 


129.8 


186.4 


2.20 


C/1993 Fl 


59 - 


1 5.9005 


-45.0 


330.2 


213.7 


1.44 


C/2002 J5 


60 - 


1 5.7268 


-33.0 


240.2 


157.1 


67.57 


C/2005 L3 


60 - 


1 5.5933 


-30.3 


203.9 


169.6 


61.07 


C/2000 Yl 


61 - 


1 7.9748 


28.0 


350.3 


8.5 


33.52 


C/1999 H3 


65 - 


1 3.5009 


-50.1 


250.1 


167.3 


49.12 


C/1898 LI 


68 - 


1 1.7016 


28.3 


143.8 


315.4 


2.01 


C/1999 N4 


68 - 


1 5.5047 


-23.3 


200.6 


204.1 


61.14 


C/1972 LI 


69 


1 4.2757 


-40.3 


235.4 


66.6 


62.09 


C/2007 JA21 


69 - 


1 5.3682 


-31.4 


272.5 


247.5 


46.07 


C/1958 Rl 


76 - 


1 1.6282 


-18.4 


307.6 


133.4 


14.74 


C/1955 Gl 


82 - 


1 4.4957 


-29.8 


233.2 


168.5 


72.73 


C/2005 Al-A 


94 - 


1 0.9069 


29.4 


113.9 


12.6 


28.06 
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